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Abstract 

Propagation of dressed solitary excitations are studied in a partially degenerate quantum plasma 
in the framework of quantum-hydrodynamics (QHD) model using multiple scales technique. The 
evolution equation together with a linear inhomogeneous differential equation is solved using 
Kodama-Taniuti renormalizing technique. It is shown that the type of solitary excitations (bright 
or dark) is defined by two critical plasma parameter values. 
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I. INTRODUCTION 



Quantum effects are of vital importance on collective phenomena in dense plasmas |l| 
which may be encountered in astrophysical objects like white dwarfs, active galactic nuclei 
and neutron stars ^ or in metals and artificially manufactured compounds like semicon- 
ductors and nano-structured materials 3|, One of the key features of dense plasmas is 
the quantum degeneracy effect due to Pauli exclusion principal which starts playing an ef- 
fective role when the de Broglie thermal wavelength \b = /i/(27rmJcB7')^^^ of the plasma 
ingredients is larger than the average inter-particle distances, n~^/^ jsj or equivalently when 
the particles temperature is below the characteristic Fermi-temperature. This condition is 
well satisfied for metals and dense astrophysical objects. On the other hand, a zero temper- 
ature Fermi gas model, in which all electrons are in degenerated state, may be an idealistic 
assumption for real hot quantum plasmas such as white dwarfs j^. Therefore, at some cir- 
cumstances such as the white dwarf surface one may expect two types of electrons, namely, 
degenerate and non-degenerate. 

On the other hand, recently the field of quantum plasmas, has achieved much attention 
due to diverse applications in rapidly growing miniaturization technology [l-lS]. Different 
aspects of nonlinear wave phenomenon in quantum plasma systems has been studied using 
quantum hydrodynamics (QHD) model 14j-|l7l|. Also, higher nonlinearity effects which is 
supposed to play significant role in nonlinear wave dynamics in quantum plasmas has been 
investigated in electron-positron-ion plasmas, where it has been shown that the quantum 



tunneling effects, caused due to quantum force, gives rise to higher dispersion rather t 



nan 



dissipation of solitary propagations and greatly affects the nonlinear wave dynamics 



18| 



The small amplitude electron-acoustic solitary excitations in partial 
;emperature) quantum plasmas has been studies by number of authors 



degenerate (two- 
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22|. Misra et.al 



23| have found ranges of cold-to-hot electron density ratio and the quantum diffraction 



parameter values which correspond to different types of solitary excitations. Also, recent 
investigation of electron-acoustic solitary propagations has shown that critical hot-to-cold 
electron densities can exist in two electron species quantum plasmas 2J]. On the other hand, 
it has been noted that the propagation of the electron-acoustic waves remain undamped in 
the range 0.25 < rico/riho < 4 [25|], where, Ucq and Uho denote the cold and hot electrons 
equilibrium densities, respectively. 
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In current study the QHD model is applied to investigate the higher-nonlinearity effects 
on electron-acoustic solitary excitations in two-electron-type quantum electron-ion plasma. 
The results obtained here cab be helpful in astrophysical dense-plasma such as white dwarfs. 
The organization of the article is as follows. Description of quantum hydrodynamics state 
equations is given in Section II. Reductive perturbation method is applied and the Korteweg 
de Vries (KdV) evolution equation is obtained in Section III. The stationary solution to 
higher-order soliton amplitude approximation is given in Section IV. Section V presents the 
discussions based on the numerical analysis and, finally, Section VI devotes to the concluding 
remarks. 



II. QUANTUM HYDRODYNAMICS MODEL 



In this section we consider a collision-less unmagnetized quantum plasma with two types 
of electrons, namely, degenerate and non-degenerate and additional neutralizing background 
inertial ions. This model may be appropriate for characterizing a plasma near the gaseous 
fringe of a white dwarf 26| , where both degenerate and non- degenerate electron populations 
may be present. The degeneracy pressure only applies to hot ions with temperature below 
the corresponding Fermi temperature. On the other hand, the collisions in such plasma are 
assumed to be limited due to Fermi-blocking process, hence, to be considered as collision-less. 
The closed set of basic quantum hydrodynamics (QHD) equations, involving the potential 
tunneling force, can be written as 
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We used the subscripts n and d to identify the non-degenerate and degenerate plasma 
species, respectively. The parameter h is the normalized Plank constant and Ni and Zi 
denote the background-ion density and charge. It is noted that, Z„ = = 1 and = m^. 
From the standard definitions the degeneracy pressure relates to the degenerate electrons 
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number-density through the following relation 



Pd = o 2 ^rf' '^'^rf = \/ ' = ksTpd, (2) 

OTi^Q V '^d 

where, quantities vfcL-, Epd and Tpd denote Fermi-velocity, Fermi-energy and Fermi- 
temperature, respectively and n^ o is the degenerate electrons equilibrium number- density. 
Under the zero-temperature assumption for degenerated electrons the Fermi-temperature is 
related to the number-density via the relation 

7^™=^ = ^(3-VoP- (3) 

On the other hand, we use the following scalings together with Eqs. ([3]) and ([2]) to obtain 
a dimensionless set of QHD equations, as 

Csd , t 2kBTpd , . 

X — X, t — > , n — )■ nudfi, u — )■ uCsd, ^ — ^ ^ • (.4) 

^pd ^pd ' e 

where, Upd = a/ e^ndfi/somd and Csd = ^/'^kBTpd/md are the characteristic quantum plas- 
mon frequency and Fermi-speed, respectively. The normalized set of equations read as 



dt dx ^' 

at dx 

dug _|_ dud — d(j>_ _ dnd _i_ H'^ d 
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a^n _|_ „, du„ — dip , d 

at ax ax ^ 2 ax 



1 a^^/n^ 



= Ud + rin- ZiNi. 

The quantities Ud (un), rid {fin) and refer to the velocity, density of degenerate (non- 
degenerate) species and the electrostatic potential, respectively. The normalized fractional 
parameter H = ^Wpd/^ksTpd is the quantum diffraction parameter, being the ratio of elec- 
tron plasmon-energy to the Fermi-energy. 

The quasi-neutrality condition is achieved via Poisson's relation at thermodynamic equi- 
librium state as 

ridfl + rinfi - ZiNi = 0, (6) 
or in dimensional form as 

1 + ^ = 6, /3 = ^, 6 = ^. (7) 
ridfi Udfi 
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The parameter /3 is the fraction of non-degenerate to degenerate electron density which can 
be considered as a measure of degree of plasma degeneracy. 

We may obtain the linear dispersion relation by a Fourier analysis of Eqs. dSJ as 

/3 , 1 



-fc^. (8) 



There are higher and lower small- limits branches in dispersion relation in the following 
forms 



_§_1.2 , 
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(9) 



t should be noted that the QHD model does not apply to large-A; values (e.g. see Ref. 



161]). The variations of dispersion curves with respect to the quantum parameter, if, and 
the degeneracy parameter, /3, is given in Fig. (1). It is clearly observed that the increase 
in the value of quantum diffraction, if, causes a larger dispersion of both higher and lower 
modes. It is also remarked that the higher mode starts at higher frequencies as the plasma 
becomes less degenerate (i.e. as /3 increases). On the other hand, at the special case of a 
classical non-degenerate plasma, /3 — )■ oo, if = the dispersion relation for acoustic mode 
reduces to ~ fc. 

III. FIRST-ORDER AMPLITUDE EVOLUTION 

Here, we introduce a strained coordinate in which the amplitude perturbation moves with 
a phase velocity of A, 

^ = e^{x- At) , (10a) 



r = eh (10b) 

Hence, the normalized QHD equations can be written in the following forms in the new 
stretched coordinate 
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with the parameter e characterized as small, positive and real number proportional to the 
perturbation amplitude. The asymptotic expansion of dependent plasma variables around 
thermodynamics equilibrium is, then, achieved using the following expansions 

nd = 1 + £n^^^ + e^n^^^ + . . . , 



Un = eu^n^ + e^U'n' + S^uf^ + . . . , 
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By isolation of the distinct perturbation orders, from the leading-orders in £, we deduce the 
following relations among first-order perturbed plasma quantities 
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Therefore, the first-order approximations for degenerate and non-degenerate ion velocities 
and densities are found to be 
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The nonlinear dispersion relation as well as the phase velocity of wave, which are obtained 
through the compatibility requirement, read as 

' f ^ 0. (15, 



A2 A 



In the next higher-order of e we get 
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Consequently, using the dispersion relation (Eq. ITSl) to eliminate the common terms in Eqs. 
( |T7|) . a KdV-type evolution equation is obtained 

%^B^a.^a^5^ = o. (18) 



The coefficients of KdV evolution equation read as 

8(l+/3f/^y^' ^ ^ 

4/3-3 

This is the equation which describes the first-order evolution of the electrostatic potential, 
however, in order to achieve better approximations we need to consider the higher terms in 
e, which will follow next. First, we note that there is critical plasma value, H^r at which the 
dispersion coefficient vanishes. Also, for the value of f3cr the nonlinear coefficient vanishes. 
Therefore, no KdV-type solitary excitations are expected for these limiting cases. To obtain 
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a solitary solution for these limiting cases one has to solve an MKdV-type equation. These 
critical plasma values are given by 



/3c 



3 
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(21) 



The higher-order perturbed components for degenerate and non-degenerate species can be 
obtained by solving Eqs. fll4p and (IT7|) together. These solutions are of the following shapes 
in terms of the first and second-order potentials for degenerate electrons 

-I- /Vioj ",'^„ . 
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and, also for non-degenerate electrons 
with the coefficients given bellow 
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Finally, by considering the third order of approximation we obtain the next-order defining 
relations for both degenerate and non-degenerate electrons as 
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"d a? "^"d '*d d(, (9? '^'''d '^v'-d ) d(, 
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Equations (!26l) together with Eqs. ( IT4l) . (|T8l) . ( l24l) and ( l25l) give rise to the following linear 
inhomogeneous equation, the solution of which along with Eq. ( IT8l) . describes the second- 
order nonlinear electron-acoustic wave amplitude evolution 
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IV. EVOLUTION OF DRESSED SOLITARY EXCITATIONS 

In order to obtain the second-order approximation to dressed amplitude evolution, the 
Eqs. ( ITSj) and (l27l) have to be solved together in a bounded form. To obtain such solution, we 
will employ the renormalizing technique introduced by Kodama and Taniuti [27]. Therefore, 
we first rewrite the Eq. (^7^ in the following form 
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where, 

r 



dr 



(30) 



Note that has be introduced for the renormahzed amphtude component. For the equation 
set describing the first and second-order amphtude evolution, we have 
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FoUowing the methodology, also described in Ref. [28|, a moving frame C = ^ ~ {M + 6M) t 
is introduced, in which Eqs. (PTj) read as 
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However, the desired bounded single-soliton stationary solutions should satisfy the following 
boundary conditions 
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We simply obtain a KdV-type solution for the first-order amplitude from Eq. fl32ap as 
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Now, using Eq. ( 130|) and the solution Eq. ( 134|) the right hand side of Eq. ( 132|) is rewritten 
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The last term in Eq. (135|) is set equal to zero to eliminate the secular term in Eq. (l32i) which 
rather simplifies Eq. fl32|) to the following form 



+ 2 (50osech2 (J) - f ) 0(2) = Cisech^ (i) - C^sech^ (J) , (36) 



where the C coefficients are given as 
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Employing the "tanh" method i.e., by changing the variable rj = tanh(|-), the Eq. ( 136|) is 
transformed into the familiar form of associated Legendre function [28]. Hence, the total 
amplitude of the perturbed electrostatic potential up to the second-order approximation, 
after absorbing the e and coefficients, can be written as 



^tot (C) = 0^'^ + 0^'^ = sech^ 



6o+('^ + ^-^sech2 
° ' 6 4 8 VA 
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V. NUMERICAL EVALUATION AND DISCUSSION 



Evaluation of the first-order amplitude solution reveals that the soliton amplitude diverges 
at the critical value of f3cr = 3/4. In this analysis we use only values of /3 ^ 1 {riQ^n ^ 
no,d) where the 0o ^ 1 is satisfied. Figure 2(a) shows two distinct regions (separated by 
the critical if-curve) where the solitary bright- or dark-type excitations can exist in this 
plasma. For the upper/lower region only bright/dark solitons exist. It is also shown in 
Fig. 2(b) that, depending on the value of matching-speed M = ±0.1, the hight/depth 
(first-order amplitude) of compressive/rarefactive solitary profile increases as the fraction of 
equilibrium nondegenerate to degenerate electron density increases. The distinct branches 
(interconnected at a value of Her) shown in Fig. 2(c) correspond to compressive (left-branch) 
and rarefactive (right-branch) solitary wave and indicate that the soliton width vanishes at 
ifcr-values. 

The first-order (dashed-curves) and the corrected (solid-curve) amplitude variations with 
respect to plasma parameters are shown in Figs. 3 and 4 for compressive and rarefactive 
solitary profiles, respectively. The left/right columns of figures show the variation with 
respect to change in /3(fixed H)/H (fixed /3). The comparison of plots indicate that for 
values of H very close to critical value (Fig. 3(d)) the correction is very large. Also it is 
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remarked that for compressive excitations the total (corrected) amphtude increases with 
increase in the value of ^ (Figs. 3(a)-3(c)), and conversely, the width of corrected profile 
tends to decrease. On the other hand, for rarefactive profiles the corrected amplitude tends 
to increase with increases of /3-values (Figs. 4(a)-4(c)), while the width tends to increase. 
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FIGURE CAPTIONS 



Figure- 1 

(Color online) Dispersion curves of linear lower and higher electron-acoustic branches 
(Eq. dH])) are shown for different values of fractional plasma degeneracy and fixed value 
of H (Fig. 1(a), 1(b)) and for different values of quantum diffraction parameter {H) and 
fixed value of /3 (Fig. 1(c), 1(d)). Different dash sizes refer to and are related appropriately 
to different values of varied parameter in each plot. 

Figure- 2 

(Color online) (a) A plot of P-H plane showing the regions where the electron-acoustic 
solitary excitations are rarefactive (lower-region) or compressive (upper-region) and (b) Vari- 
ation of the first-order soliton amplitude with respect to increase in the value of nondegen- 
erate to degenerate electron equilibrium density for fixed other plasma parameters, (c) the 
variation of soliton width with respect to quantum diffraction parameter, H showing critical 
values at which the width vanishes. Different dash sizes refer to and are related appropri- 
ately to different values of /3 in Fig. 2(b) and the left-/right-branched in this plot represent 
where the solitons are compressive/rarefactive. 

Figure- 3 

(Color online) Profiles for bright-type solitary excitations is given for first-order ampli- 
tude (dashed curves) and the corrected (total) amplitude (solid-curves) perturbations. The 
left/right columns of figures show the variation with respect to change in /3(fixed if)/if (fixed 
(5). The value of M = — 1 is used for all plots in this figure. 

Figure-4 

(Color online) Profiles for dark-type solitary excitations is given for first-order ampli- 
tude (dashed curves) and the corrected (total) amplitude (solid-curves) perturbations. The 
left/right columns of figures show the variation with respect to change in /3 (fixed if) /if (fixed 
(3). The value of M = — 1 is used for all plots in this figure. 
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0.000 
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j8=0.01,H=0.1 



yS=0.05,H=0.1 



0.000 



0.005 



0.010 



0.015 



0.020 




0.00 



-0.01 
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-0.03 



-0.04 



0.00 
-0.01 
-0.02 
-0.03 
-0.04 
-0.05 
-0.06 
-0.07 



10 



j8=0.05,H=0.1 




^=0.1,H=0.1 






(c) 



10 



0.00 



-0.01 







-0.02 



-0.03 



-0.04 








0.00 



-0.01 



-0.02 



-0.03 - 



-0.04- 



-10 







0.00 
-0.01 
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-0.03 
-0.04 
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(d) 

P=0.05, H=0.2 




yS=0.05, H=0.3 




-10 -5 
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